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ABSTRACT 

We study the multi-dimensional geometry of supernova (SN) explosions by means of spectropolari- 
metric observations of stripped-envelope SNe, i.e., SNe without a H-rich layer. We perform spectropo- 
larimetric observations of 2 stripped-envelope SNe, the Type lb SN 2009jf and the Type Ic SN 2009mi. 
Both objects show non-zero polarization at the wavelength of the strong lines. They also show a loop 
in the Stokes Q — U diagram, which indicates a non-axisymmetric, three-dimensional ion distribution 
in the ejecta. We show that five out of six stripped-envelope SNe which have been observed spectropo- 
larimetrically so far show such a loop. This implies that a three-dimensional geometry is common 
in stripped-envelope SNe. We find that stronger lines tend to show higher polarization. This effect 
is not related to the geometry, and must be corrected to compare the polarization of different lines 
or different objects. Even after the correction, however, there remains a dispersion of polarization 
degree among different objects. Such a dispersion might be caused by three-dimensional clumpy ion 
distributions viewed from different directions. 

Subject headings: supcrnovae: general — supernovae: individual (SNe 2009jf, 2009mi) — techniques: 
polarimetric 



1. INTRODUCTION 

1.1. Explosion Geometry of Supernovae 

Core-collapse supernovae (SNe) are the explosions of 
massive stars at the end of their lives. Besides their im- 
portance as the end points of stellar evolution, SNe play 
an important role on chemical and dynamical evolution 
of galaxies. However, the explosion mechanism of SNe 
ha s been unclear for a lo ng time after the first concept 
bv lBur bidae et al. (195^) and th e first numerical simula- 
tion bv lColgate fc White! (fl966h . 

Modern detailed numerical simulations agree 
that a successful explosion cannot be obtai ned 
one-dimensional s im ulatio ns dRampp &: Ja nka 
Liebendorfer et all ]2001b iThompson et al 
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attention is being pai d to multi-dimensional effects, such 
as convection (e.g.. iHerant et al.l Il994t iBurrows et al.1 
Il995t Uanka fc Muellerl I1996D and Standing Accretion 



Shock Instability (SASI, e.g., 
Scheck et all 120041: lOhnishi et all 



Blondin et all 120031: 



20061: FoglizzolEal] 



2007T llwakami et al. 1 120081) . In fact, some successful 
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explosions h aye been obtained i n two-dimensiona l (2D 
simulations dBuras et al.l 120061: iMarek fc Jankal 1200 
iSuwa et al.ll2010n . 

Thanks to the advance of numerical simulations, 
some simulations have been done in fully threc- 
dimensional (3D) s p ace 
Wongwathanarat et al.l 



Kuroda et al 



12012! ) . 
that it 



dNordhaus et al 
120101: iTakiwaki et al.l 



Interestingly, 

(|2010t ) show that it is easier to obtain a successful 
explosio n in 3D simulation than in 2D simulations (but 
see also lHanke et aLlfeOllh . 

Given these circumstances, it is important to obtain 
observational constraints on the multi-dimensional ge- 
ometry of SN explosions. The most direct way is the ob- 
servat ion o f spatially resolved su pernovae. Ilsensee et al.1 
(|2010l ) and lDeLanev et aD (|2010[ ) performed detailed ob- 
servations of Galactic supernova remnant Cassiopeia A. 
Using a combination of imaging and spatially resolved 
spectroscopy, they reconstruc ted the 3D geome try of the 
SN. Similarly, direct imaging (IWang et alj|2002D and spa- 
tially resolved spectroscopy (|Kiaer et al.ll20101 V have been 
performed for SN 1987A in the Large Magellanic Cloud. 
However, the number of objects suitable for such obser- 
vations is limited. 

Observations of extragalactic SNe can also provide 
hints of explosion geometry, although at this dis- 
tance SNe are point sources. For example, spectro- 
scopic observations of SNe at >1 yr after the explo- 
sion provide information of explosion geometry. At 
such late epochs, SN ejecta are optically thin. There- 
fore, this method is sensitive to the geometry of the 
innermost ejecta. By analyzing line profiles of nebu- 
lar emission, bipolar geometry of SNe has been sug- 
gested (e.g., iMazzali et all 120011 IMaeda et all l2002t 

ot oi Ion 
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.211121: 
iNordhaus et al.1 
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Mazzali et al.l 120051: IMaeda et al.1 120081; iModiaz et all 
20081: iTanaka et all l2009bt 



Tauben berger et al.1 12009^ 
Maurer et al.l2010l but see also lMilisavlievic et al.l2010l ). 
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(b) Spherical ion distribution 



(a) Spherical photosphere 




Note that an asymmetric ion distribution is ex- 
pected either from an asymmetric element distribution, 
which could be produced by a multi-dimensional explo- 
sion (e.g., iNaeataki et all Il997t iKifonidis et~ai1 120001 



(c) Aspherical ion distribution 



I Polarization vector 
(length = degree of polarization) 

Fig. 1. — S chematic illustration of polarization in the SN 
ejecta (see also I Kasen et al . 2003; Leonard k, Filippenko [20051 ; 
Wang & Wheeler 2008^ (a) When the photosphere is spherical, 
polarization is cancelled out, and no polarization is expected. At 
the wavelength of a line, polarization produced by the electron 
scattering is depolarized by the line transition, (b) When the ion 
distribution is spherical, the remaining polarization is cancelled, 
and no polarization is expected, (c) When the ion distribution is 
not spherical, the cancellation be comes incomplete, and line po- 
larization could be dete cted (e.g., IJef fcrv 1989; IKasen et al. 2003; 
IWang fc W heeler 2003). 



Polarization at early phases (< 50 d after the ex- 
plosion) is another powerful method to study the 
multi-dimensional geom etry of extragalactic SNe (see 
IWang fc WlieeTei1[2008l for a review) . In contrast to late 
phase spectroscopy, polarization is sensitive to the geom- 
etry of the outer ejecta. In particular, spectropolarime- 
try can probe 3D geometry, as explained in the following 
section. In this paper, we study the multi-dimensional 
explosion geometry of SNe using spectropolarimetric ob- 
servations. 

1.2. Power of Spectropolarimetry 

Polarization in SNe is generated by electron scatter- 
ing. Since the polarization degree of the scattered light 
becomes maximum when the scattering angle is 90°, the 
expected polarization map for a spherical photosphere 
looks like (a) in Figure [T] However, SNe are observed as 
point sources, and thus all the polarization vectors cancel 
out, resulting in no polarization. 10 

Absorption and re-emission by a bound-bound tran- 
sition, which make a P-Cygni profile in the expand- 
ing ejecta, tend t o cause depolarization (jJeffervi fl989r , 
IKasen et al.l l2003h . Thus, the polarization produced by 
the electron scattering is effectively concealed by the in- 
teraction with lines. If an ion is distributed spherically 
(b), polarization vectors are still cancelled out, and no 
polarization is expected at the wavelength of the line. 
Polarization would be detected at the wavelength of the 
absorption line only if the ion is distributed asymmetri- 
cally (c), because of the incomplete cancellation. Thus, 
we can identify an asymmetric ion distribution by de- 
tecting line polarization. 

10 When the photosphere is deformed, cancellation of polariza- 
tion become s incomplete. As a result, continuum polar ization will 
be detected ^Shapiro fc Sutherland|[l98 2; Hoflich 1991;!). Although 
continuum polarization is difficult to distinguish from interstellar 
polarization, multi-epoch observations have revealed that there is 
non-ze ro continuum polarization in Type UP SNe after the plateau 
phase ULeonard et alj|2006t IChornock et aI.II2010r) . 



20031: iMaeda et al.ll2003t lTominagall2009l : iFuiimoto et al.l 
201 lh . or by an asymmetric ionization/temperature 
structur e, which is caused b y an asymmetric distribution 
of 56 Ni (|Tanaka et al.l [2001 . 

The properties of linear polarization are fully described 
by the Stokes parameters. Throughout this paper, we de- 
fine the Stokes parameters as a fraction of the total flux: 
Q = Q/I and U = U/I, where Q and U can be expressed 
by Q = Iq — /90 and U = I45 — -Z135, respectively. Here Ie 
is the intensity measured through the ideal polarization 
filter with an angle 9. From the Stokes parameters Q 
and U, the polarization position angle 9 is obtained by 



29 = atan(C//Q). 



(1) 



The position angle is measured from north to east on the 
sky. 

Figure [2] shows the expected line polarization in a 2D 
bipolar model (top) and a 3D clumpy model (bottom). 
The polarization spectra are calculated by Monte Carlo 
method for arbitrary line optical de pth distributions, un - 
der assumptions similar to those in IKasen et al.l (|2003h : 
iHole et aLT (|2010[ ) . The details of the simulations will be 
given in a forthcoming paper. The left panels show the 
distribution of line optical depth. In the 2D model, the 
opacity is enhanced by a factor of 10 in the polar region, 
while in the 3D model, the opacity is similarly enhanced 
in the randomly distributed clumps. 

In both cases, the polarization changes at the wave- 
length of the line (middle panels) , which is caused by the 
asymmetric distribution of line optical depth. However, 
the behavior of the polarization angle is different. In the 
2D case, the polarization angle is constant (some fluctu- 
ation is caused by the Monte Carlo noise), while in the 
3D case, it largely changes across the line. Since the SN 
ejecta expand homologously and the velocity can be used 
as a radial coordinate, a change in the polarization angle 
means a radially-dependent ion distribution. 

This effect is more clearly seen in the Q — U plane. 
The right panels in Figure [2] show the same simulated 
polarization spectrum in the Q — U plane. Different col- 
ors represent different Doppler velocity, according to the 
color bar above the plots. The polarization for the 2D 
model shows a straight line in the Q — U plane, while that 
for the 3D model shows a "loop" , which corresponds to 
a change in the polarization angle. 

Therefore, a change in the polarization angle or the 
loop in the Q — U plane is indicative of a 3D ion 
distribution. In order to fully utilize the power of 
spectropolarimetry, we should obtain high-quality data 
to discriminate such a feature from noise. Loops in 
the Q — U pl ane have been commonly observed in 
Type la SNe (IWang et~aTl l2Q03a|; IKasen et al] 120031 : 
IChornock fc Filippenko! 120081: iPatat et al.l l2009h and in 
some core-collapse SNe (e.g., iMaund et al.1 l2007bl E 
120091) . 

1.3. This Paper 

In this paper, we show our spectropolarimetric ob- 
servations of 2 Type Ib/c SNe with Subaru telescope 



3D Explosion Geometry of Stripped-Envelope Core-Collapse SNe I 



3 




y (1000 km/s) 




-25 -20 -15 -10 -5 5 

Doppler velocity (10 3 km s~ 1 ) 



-20 -15 -10 
Doppler Velocity (10 3 



-5 

km s _1 ) 



& 





Q (%) 




y (1000 km/s) 




-25 -20 -15 -10 -5 5 

Doppler velocity (10 3 km s~ 1 



Fig. 2. — Simulated line polarization for a 2D bipolar model (top) and 3D clumpy model (bottom). (Left) The distribution of the line 
optical depth. (Middle) The total flux spectrum and polarization spectrum. The polarization angle is consistent with a constant in the 2D 
model, while it changes across the line in the 3D model. (Right) The same simulated polarization but in the Q — U diagram. Different colors 
represent different Doppler velocity, according to the color bar above the plots. The polarization data show a straight line (constant angle) 
in the 2D model, while the data show a loop (variable angle) in the 3D model. Details of the simulations will be given in a forthcoming 
paper. 



(Sections Hand [3]). The first targ et is the Type lb SN 
2009jf (discovered in NGC 7479 bv lLi et al.ll2009l on UT 
2009 September 27.33, and cla ssified by Kasliwal et all 
2001 iSahu et al.1 f2009h . See iSahu et al.1 (|201lD and 
Valenti et al. I (1201 if for more details. The other tar- 
get is_Jhe_^YMlc_SN 2009mi (discovered in IC 2151 
bv lMonard et al.ll2009l on UT 200 9 December 12.91, and 
classified bv lKinugasa et ail l2009t ) . 

We then summarize the spectropolarimetric data of 
stripped-envelope SNe published so far (Section 0}. We 
argue that a non-axisymmetric, 3D geometry is common 
in stripped-envelope SNe. We find a relation between line 
polarization and the depth of absorptions, i.e., a stronger 
line tends to show a higher line polarization. Even after 
correcting for the effect of the absorption depth, there 
remains a dispersion in the polarization degree among 
different objects. The implications of this dispersion are 
discussed. Finally, we give conclusions in Section [S] 

2. OBSERVATIONS 

We have performed spectropolarimetric observations 
of SNe 2009jf and 2009mi with the 8.2m Subaru tele- 



scope equipped wit h the Faint Object Cam era and Spec- 
trograph (FOCAS. iKashikawa et all[200l on UT 2009 
October 24.3 (MJD=55128.3) and 2010 January 8.3 
(MJD=55204.3), respectively. These epochs correspond 
to t = +9.3 and +26.5 days from the B ban d maximum 
(MJD =55118.96 for SN 2009jf according to ISahu et afl 
[20lH and MJD = 55177.8 for SN 2009mi, based on our 
observations). Hereafter, t denotes the days after the B- 
band maximum. The log of observations are shown in 
Table HJ 

For both observations, we used an offset slit of 0.8" 
width, a 300 lines mm -1 grism, and the Y47 filter. This 
configuration gives a wavelength coverage of 4700-9000 
A. The wavelength resolution is AA ~ 10 A. For the 
measurement of linear polarization, we use a rotating su- 
perachromatic half-wave plate and a crystal quartz Wol- 
laston prism. One set of observations consists of the in- 
tegrations with 0°, 45°, 22.5°, and 67.5° positions of 
the half- wave plate. From this one set of exposures, 
Stokes parameter s Q and U are derived as described by 
iTinbergenl (|1996h . 

For the observations of SN 2009jf we performed six 
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TABLE 1 

Log of Observations 



Object 


Date (UT) 


Date (MJD) 


Exposure time (s) 


Airmass 


Comment 


SN 2009jf (+9.3 d) 


2009 Oct 24.3 


55128.3 


(600 X 4) X 6 


1.03 - 1.43 


SN 


BD+28°4211 


2009 Oct 24.2 


55128.2 


(20 X 4) + (40 X 4) 


1.01 


unpolarizcd/flux std. 


G191-B2B 


2009 Oct 24.6 


55128.6 


(60 x 4) x 2 


1.33 


unpolarizcd std. 


Hiltncr 960 


2009 Oct 24.2 


55128.2 


20 x 4 


1.08 


polarized std. 


SN 2009mi (+26.5 d) 


2010 Jan 8.3 


55204.3 


(600 x 4) + (1000 x 4) 


1.26 - 1.83 


SN 


G191-B2B 


2010 Jan 8.2 


55204.2 


(60 x 4) x 2 


1.43 


unpolarizcd/flux std. 


HD 14069 


2010 Jan 8.2 


55204.2 


(5 x 4) x 2 


1.03 


unpolarizcd std. 


HD 251204 


2010 Jan 8.2 


55204.2 


20 x 4 


1.48 


polarized std. 



Note. — All the observations were performed with a 0.8" width offset slit, a 300 lines mm 1 grism, and the Y47 
filter, giving the wavelength coverage 4700-9000 A and the wavelength resolution AA ~ 10 A. 



sets of the integrations at the four angles with a total 
exposure time 4.0 hr. For SN 2009mi, we performed five 
sets of the integrations at the four angles with a total 
exposure time of 3.8 hr. Typical seeing during the ob- 
servations was 0.8" and 1.1 - 1.5" for SNe 2009jf and 
2009mi, respectively. 

After deriving the Stokes parameters for each set of ob- 
servations, these were combined. Then, the instrumental 
polarization (~ 0-0.4 %) was evaluated and subtracted 
using unpolarized stars (| Schmidt et al.|[l992T) ; G191-B2B 
(for both objects), BD+28°4211 (for SN 2009jf) and 
HD 14069 (for SN 2009mi) observed on the same night. 
The reference axis of the position angle was calibrated 
by the o bservation of the str ongly polarized star Hilt- 
ner 960 dSchmidt et all Il992l for SN 2009jf) and HD 
251204 ( jTurnshek et alJl990T for SN 2009mi). The wave- 
length dependence of the optical axis of the half-wave 
plate was corrected using the dome flat taken through a 
fully-polarizing hlter. The total flux was calibrated us- 
ing the observation of the spectrop hotometric standard 
stars BD+28°4211 and G191-B2B jMassev etall 119881 : 
iMassev fc GronwaUlToM fctk^fTOQCTh for SNe 2009jf and 
2009mi, respectively. Telluric absorption lines were also 
removed using the spectrum of these standard stars. 

3. RESULTS 

Figure [3] shows the observed total flux and polariza- 
tion spectra of SNe 2009jf and 2009mi. The observed 
polarization is the sum of the intrinsic and interstel- 
lar polarization (ISP). ISP is caused by interstellar dust 
both in our Galaxy and in the host galaxy. However, 
since the wavelength d ependence of the ISP is smooth 
(Scrkows ki et al.lll975| ). the changes in the polarization 
at the strong absorption lines, such as He I, O I, and 
Ca II, are certainly intrinsic. These features clearly in- 
dicate a broken symmetry of the SNe. 

The same data are shown in the Q — U diagram in Fig- 
ure SJ Different colors represent different wavelengths, 
according to the color bar above the plots. Most of 
the data points, except for the strong lines, are located 
around (Q , U ) ~ (-0.25 %, 0.50 %) and (-0.30 %, 
0.60 %) for SNe 2009jf and 2009mi, respectively. Spec- 
tropolari metric type is Nl accordi ng to the classification 
scheme of lWang fc Wheeler! (|2008| ) . Some deviation from 
the reference point could be due to non-zero continuum 
polarization of SNe. Hereafter we use these (Qo, Uo) as a 
reference point. This component is most likely to be the 
ISP for the following reason. The ISP can be corrected if 
we assume a complete depolarization at strong emission 



lines (e.g., 


Kawabata et al.l l2002t iLeonard et al.l 120051: 


IWang et al.l 


20061: iMaund et al.l2007b.a). This method 



is justified because the emission part of a P Cygni profile 
consists largely of line-scattered, depolarized photons. In 
our data, the polarization at the strong emission peak at 
the Ca II line is consistent with the reference point above. 
In the following discussion, however, we do not explicitly 
assume this component to be the ISP but focus only on 
the line polarization. 

The Q — U diagram (Figure [4j) is a more unambigu- 
ous way to show the observed polarization data than a 
plot as a function of wavelength (Figure ■ Since the 
ISP could add an offset to the SN polarization, the po- 
larization angle 6 in Figure [3] is not directly related to 
the SN properties. Instead, one can readily measure the 
polarization angle 9' from the reference point (Qo, Uq); 

2<?' = atan[([/-[/ )/(Q-Q )]. (2) 

In Figure [5j we show the data only around the O I 
and Ca II lines as a function of Dopplcr velocity. It is 
clear that the data at the Ca II and O I lines occupy 
different regions in the Q — U diagram, indicating differ- 
ent spatial distributions between Ca II and O I. Such a 
difference i s also clearly seen i n other SNe e.g., Type lb 
SN 2008D (IMaund et all [2001) and Type lib SN 2008ax 
(|Chornock et al.ll2011lh 

A more interesting feature is the shape of the polar- 
ization data in the Q — U diagram. Starting from the 
reference point, the Ca II and O I lines in SN 2009jf 
show a loop at these lines (spectropolarimetric type L). 
This means that the angle 9' measured from the refer- 
ence varies with Doppler velocity, and the depth in the 
ejecta (homologous expansion, r = vt). As shown i n 
Figured and as sug gested by e.g., iKasen et al.l (|2003l) : 
IMaund et al.l (|2007bl fq). this loop clearly indicates that 
even axisymmetry is broken in the SN ejecta. 

The Ca II feature of SN 2009mi is even more intriguing. 
Measuring from the reference point, the variation in the 
angle 9' is as large as ~ 90° (the difference of 180° in the 
Q — U diagram corresponds to the difference of 90° in 
the polarization angle on the sky) . This indicates a large 
change in the distribution of the Ca II ion depending on 
the depth of the ejecta. 

4. PROPERTIES OF SUPERNOVA LINE 
POLARIZATION 

4.1. Ubiquity of Non- axisymmetry 

In Table [H we summarize line polarization measure- 
ments obtained so far for stripped-envelope SNe. Here 
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Fig. 3. — Total flux spectrum and polarization spectrum of SNe 2009jf (left) and SN 2009mi (right). The total flux is shown in units 
of 10~ 15 and 10" 16 erg s 1 cm 2 A for SNe 2009jf and 2009mi, respectively. The polarization data are binned into 25 and 50 A for SNe 
2009jf and 2009mi, respectively. In the plot of Stokes parameters and polarization angle, contribution of the interstellar polarization is not 
corrected for. 
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Fig. 4. — Observed polarization data of SNe 2009jf (left) and 2009mi (right) in the Q — U plane. The polarization data are binned 
into 50 A. The interstellar polarization is not corrected for. Different colors show different wavelengths as shown in the color bar above the 
plots. The polarization data at the O I and Ca II lines are also shown as a function of Dopplcr velocity in Figure [5] 



the line polarization is denned as a relative change 
from the continuum level, and thus, not affected by the 
ISP. We show 6 stripped-envelope S Ne with high-quality 
data, namely Type lb SNe 2 005bf (IMaund et al l l2007bl: 
iTanaka et aLll2009al) . 2008D (|Maund et al.ll2009fl. 2009if 
(this paper) , Type Ic SNe 2002ap dKawabata et al.l 
[200l iLeonard et al l [20021: IWang et al.1 l2003bD . 2007ct 
(ITanaka et al.ll2008l), and 2 009m i (this paper). As also 
noted bv IWang fc Wheeler] (|2008f) . it is clear that all SNe 
show non-zero polarization, which means that stripped- 



envelope SNe generally have asymmetric explosion geom- 
etry. 

For SNe 2005bf and 2008D. IMaund et al.l (|2007Ul2009l) 
found a loop in the Q — U diagra m at strong lines , 
which is indicative o f a 3D geometry ([Kasen et al.ll2003l : 
IMaund et all I2007bllcl . see also Figure [5]). Our new 
data for SNe 2009jf and 2009mi also show a loop in 
the Q — U diagra m. We have checked the literature 
abou t SN 2002ap ifKawabata et al.ll2002l ; [Leonard et al.l 
I2002t IWang et all l2003bl) . Although they do not ex- 
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Fig. 5. — Observed polarization data of SNe 2009jf (left) and 2009mi (right) around the O I line (circles connected with the blue line) 
and Ca II line (squares connected with the red line). The polarization data are binned into 50 A, giving a velocity resolution of about 
1900 and 1700 km s _1 for the O I and Ca II lines, respectively. The interstellar polarization is not corrected for. Different colors show 
different Dopplcr velocities as shown in the color bar above the plots. Although the Ca II triplet seems to be resolved in the total flux of 
SNe 2009jf and 2009mi (Figure the velocity in this Figure is simply measured from the mean wavelength (8 567 A). In our data , there 
is no evidence for an additional high-velocity component of the Ca II line as seen in some other SNe (see e.g. JMaund et all 120091) . The 
O I and Ca II features show a loop in the Q — U plane, indicating non-axisymmctric distribution. In addition, the Ca II line in SN 2009mi 
shows a large change of the angle measured from the reference point. 



plicitl y mention the loop, the data in iKawabata et all 
(|2002D show a loop in the Ca n line. Thus, loops 
are quite common in stripped-envelope SNe; five of six 
stripped-envelope SNe show the loop. This implies that a 
non-axisymmetric, 3D geometry is common in stripped- 
envelope SNe. 

4.2. Relation Between Line Polarization and 
Absorption Depth 

Figure [5] shows the observed line polarization as a func- 
tion of epoch from maximum brightness. There is a large 
variety in the polarization degree. It is clear that the 
Ca II line (blue) tends to be more polarized than the 
other lines. There is no clear trend in the time evolu- 
tion; the line polarization could increase or decrease with 
time. 

Figure[7]shows the same polarization data as a function 
of the fractional depth of the absorption. The fractional 
depth (FD) is defined as 



FD = cont abs , (3) 

where i^ont and -Fkbs are the flux at the continuum and 
absorption minimum, respectively. 

Figure [7] shows that stronger lines tend to show higher 
polarization degree. In fact, there are no observational 
points at the top left of the plot, i.e., FD < 0.3 and 
-Fobs > 1%. This seems natural because the absorption, 
unless it is quite strong, cannot cause a significant in- 
complete cancellation of the polarization (Figure [T]). 

We can analyze this relation in a simple, analytic way. 
We assume a spherical photosphere projecting an area S 
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Fig. 6. — Observed line polarization as a function of SN epoch. 
Data for the same SN are connected with lines. There is no clear 
trend in the time evolution of the line polarization. The Ca II line 
tends to be more polarized than other lines. 



on the sky. Within this photospheric disk, the flux per 
unit area / and the polarization degree are assumed to 
be uniform. The total flux from the photospheric disk is 
-Fcont = IS. When light is emitted from a point on the 
photospheric disk on the sky, it can undergo line absorp- 
tion with an absorption fraction a; a b s = 1 — cxp(— r). At 
the wavelength of this line, the flux is -Fkbs = (1— x&bs)IS- 
If the distribution of the line opacity is uniform, no 
polarization is expected at the line, i.e., P aDS = (Section 
II. 2p . To introduce asymmetry in the line opacity, we 
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Fig. 7. — Observed line polarization as a function of fractional 
depth of the absorption. Data for the same SN are connected with 
lines. There is a general trend that stronger lines show higher 
polarization, as expected from Equation [8] The gray dashed line 
shows the line polarization as a function of FD which is expected 
by Equation [8] for the case of P cor = 2.0%. 



assume that the absorption fraction is enhanced by a 
factor of / in a region AS (e.g., the shaded region in 
Figure 1(c)). Then, the total flux at the wavelength of 
the line is 

Fabs = (1 - x,hs)I{S - AS) + (1 - fx ahs )IAS 

= [(1 - ^abs) - (/ - l)x ahs AS/S]IS, (4) 

and the fractional depth (Equation [3]) can be written by 
using Equation 0] 



FD = x abs + (/ - l) Xa . hs AS/S. 



(5) 



The asymmetry introduced above results in an incom- 
plete cancellation of the polarization. If this region has 
a constant polarization direction, the amount of non- 
cancelled flux is equivalent to the excess absorption in 
the region, i.e., 



■Pabs-Fabs = (/ - l^abs/AS. 



(6) 



Then, the polarization degree can be expressed as a func- 
tion of the fractional depth by Equations [5] and [6] 



P 



abs 



if - l) Xllhs AS/S 
1 - FD 



(7) 



Unless the enhanced opacity dominates the absorption, 
FD ~ x a bs (Equation [5]) . Thus, 

FD FD 

P abs ~ (/ - 1)AS/S T -^ = ^corT^^, (8) 

where we define a corrected polarization degree P cor = 
(f—l)AS/S. This corrected polarization degree is equiv- 
alent to the polarization with FD = 0.5. It can be inter- 
preted as a rough indicator of the product of the amount 
of the enhancement (/ — 1) and the fractional area of 
the enhanced region AS/S under the assumption that 
the enhanced region only has one polarization direction. 
The gray dashed line in Figure [7] shows the line polariza- 
tion expected by Equation [5] with P cor = 2.0%. 

Note that the analytic equations here involve signifi- 
cant simplifications, i.e., a constant intensity in the pho- 
tosphcric disk, and a constant direction of polarization in 



the enhanced region. But if we use AS as an intensity- 
weighted area, in which polarization is not cancelled, the 
equation holds even for general cases with a non-constant 
intensity and many opacity-enhanced regions. 

From Equation[8j it is expected that stronger lines with 
a larger fractional depth tend to have higher polariza- 
tion, which is actually found in Figure [7J It should be 
emphasized that the variety of the observed line polariza- 
tion is largely generated by a variety of the absorption 
depth. The absorption depth is primarily determined 
by the density, chemical abundance, or ionization in the 
ejecta, which can be varied even in the ID explosion (re- 
flecting the mass and kinetic energy of the ejecta, and 
the temperature in the ejecta). 

Figure [8] shows the observed polarization (left) and po- 
larization corrected for the absorption depth (right) for 
each stripped-envelope SN. It is obvious that after correc- 
tion for the absorption depth, the Ca II line is not always 
more polarized than the other lines. This demonstrates 
the importance of the correction. To compare the polar- 
ization of different lines or different objects correctly, one 
must correct for this effect, which is not directly related 
to the geometry. This fact has been overlooked in previ- 
ous analyses. It would be interesting to study this effect 
also on the line polarization of Type la SNe, which has 
been suggested to be correlated wit h the declining rate 
of the light curve ()Wang et al.ll2007l) . 

4.3. Properties of Corrected Polarization 

It is clear that there is a remaining dispersion in the 
line polarization among different objects even after the 
correction for the absorption depth (right panel in Fig- 
ure [8]). We first examine if the corrected polarization is 
related with the epoch. The left panel of Figure shows 
the corrected line polarization function of epoch. 
Although there is a weak trend that the polarization de- 
creases with time, it is not clear if all the variety in the 
corrected polarization can be explained by the effect of 
the epoch. 

The remaining dispersion could also be understood 
by the multi-dimensional properties of each SN. It 
is thus worth studying a relation between the line 
polarization and the [O i] line profile in the neb- 
ular spectra, which is another geometric probe of 
the SN ejecta (iMazzali et al J l200¥ IMaeda et al.1 120021: 
Mazzali etaTI 120051: IMaeda et al.l 120081: iModiaz et ail 
20081: ITanaka et all l2009bl iTaubenberger et al.l 120091: 
Maurer et al.ll2010D . 



The right panel of Figure [§] shows the relation be- 
tween corrected line polarization and the [O i] line pro- 
file. SNe 2002ap and 2007gr show a single-peak pro- 
file ([Leonard et al.ll2002l: iFolev et al.ll2003t iMazzali et al.1 
I2007t iHunter et al.l l2009t iMazzali et al.l I2010D while 
SNe 2005bf and 2 08D show a doub l e-peak profile 
(IMaeda et al.l l2007t iModjaz et all 120091 : ITanaka et all 
I2009H) . The[Q i] profile of SN 2009jf is compl ex; it shows 
both single- and double -peak components ()Sahu et al.l 
201lUValenti et al.ll20Tl . Following the classification by 
Maeda et al.1 (|2008| ). we classify it as "transition". Only 



with five objects, it is difficult to draw a firm conclusion. 
It is therefore important to increase spectropolarimetric 
samples to further study this relation. 

The dispersion may also reflect a difference between 
Type lb and Ic. In fact, 3 Type Ic SNe in our sample 
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Fig. 8. — (Left) Observed line polarization for stripped-envelope SNe. The Ca II line is usually more polarized than the other lines. 
(Right) Line polarization corrected for the absorption depth (Equation [8} . After correcting the effect of the absorption depth, the Ca II 
lines is not always more polarized than the other lines. 
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Fig. 9. — (Left) Corrected line polarization as a function of SN epoch. Data for the same SN are connected with lines. (Bight) Relation 
between corrected line polarization and nebular [O i] line profile. 



tend to have lower polarization degree than 3 Type lb 
SNe (left panel of Figure [8]). By increasing the number 
of samples, we will be able to study the effect of the 
existence of the He envelope on the explosion geometry 
One natural scenario is that the dispersion is caused 
by random line-of-sight effects. As discussed in Section 
[3l stripped-envelope SNe generally have 3D geometry 
If we assume a clumpy structure as in the left bottom 
panel in Figure O a variety of polarization is naturally 
expected for different linc-of-sight. Since the effect of the 
line-of-sight is very sensitive to the size and the number 
of clumps, it can be used as an observational probe to 
study a typical 3D geometry of SNe. This will be further 
studied in a forthcoming paper. 

5. CONCLUSIONS 

We have performed spectropolarimetric observations of 
2 stripped-envelope SNe, Type lb SN 2009jf and Type Ic 
SN 2009mi. Both objects show non-zero line polarization 
at the wavelength of strong absorption lines. The polar- 
ization data show a " loop" in the Q — U diagram, which 



indicates a non-axisymmetric, 3D distribution of ions in 
the SN ejecta. After adding our new data to the sample 
of stripped-envelope SNe with high-quality spectropolari- 
metric observations, five SNe out of six show a loop in the 
Q — U diagram. We conclude that a non-axisymmetric, 
3D geometry is common in stripped-envelope SNe. 

We have studied the properties of line polarization in 
stripped-envelope SNe. We have found that a stronger 
line tends to show a higher polarization. We give an 
analytic equation that naturally explains this relation. 
Although the Ca II line usually shows a higher line po- 
larization than the other lines, it is not always true after 
the correction of the absorption depth. This effect must 
be corrected in order to compare the polarization of dif- 
ferent lines or different objects. 

Even after the correction for the absorption depth, 
there is a dispersion in the polarization degree among 
different objects. This dispersion could simply reflect the 
epoch of the observations. It might also be related to the 
[O i] line profile in the nebular spectra or SN types (Type 
lb or Ic) , although it is not possible to draw any firm con- 
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elusion with the current sample. Also, if we assume a 3D 
clumpy structure, the variety is naturally expected by 
the line-of-sight effect. Modeling of spectropolarimetric 
data will be performed in a forthcoming paper. 
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TABLE 2 
Summary of Line Polarization 
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Type 


3D? 


Epoch 


^PcII 




^OI 


-PNal/HoI 


FD FeII 


FD CaII 


FDqi 


FD NaI / HeI 


Ref. 








(day) 


(%) 


(%) 


(%) 


(%) 










SN 2005bf 


lb 


VPS 


-6 


0.8 (0.2) 


3.5 (0.5) 


0.0 (0.3) 


1.2 (0.2) 


0.37 (0.03) 


0.50 (0.05) 


0.0 (0.0) 


0.40 (0.03) 


1 
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0.4 (0.2) 


1.5 (0.3) 


0.0 (0.5) 


0.6 (0.4) 


0.21 (0.02) 


0.39 (0.02) 


0.0 (0.0) 


0.46 (0.02) 


2 


SN 2008D 


lb 


yes 


3.3 


0.8 (0.2) 


1.8 (0.3) 


0.5 (0.13) 


0.4 (0.2) 


0.35 (0.03) 


0.49 (0.02) 


0.23(0.03) 


0.46 (0.02) 


3 








18.3 


1.0 (0.3) 


2.5 (0.7) 


0.3 (0.13) 


1.1 (0.1) 


0.50 (0.05) 


0.58 (0.05) 


0.24 (0.02) 


0.58 (0.02) 


3 


SN 2009jf 


lb 


yes 


9.3 


0.4 (0.2) 


1.2 (0.2) 


0.9 (0.2) 


0.5 (0.2) 


0.32 (0.03) 


0.69 (0.02) 


0.38 (0.02) 


0.38 (0.02) 


this paper 


SN 2002ap 


Ic 


yes 


1 


0.18 (0.05) 




0.8 (0.1) 


0.0 (0.05) 


0.20 (0.03) 




0.38 (0.02) 


0.04 (0.02) 


4,5 
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0.12 (0.05) 




0.6 (0.1) 


0.0 (0.1) 


0.21 (0.03) 




0.34 (0.02) 


0.04 (0.02) 


1.5 
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0.0 (0.1) 


0.3 (0.1) 


0.5 (0.1) 


0.0 (0.1) 


0.27 (0.05) 


0.36 (0.05) 


0.50 (0.05) 


0.06 (0.02) 
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27 




1.6 (0.1) 


0.0 (0.2) 


0.2 (0.1) 




0.65 (0.04) 


0.41 (0.05) 


0.31 (0.03) 


4 


SN 2007gr 


Ic 


no 


21 


0.0 (0.3) 


2.5 (0.3) 


0.0 (0.3) 


0.0 (0.3) 


0.42 (0.03) 


0.88 (0.02) 


0.56 (0.02) 


0.65 (0.02) 
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SN 2009mi 


Ic 


yes 


26.5 


0.5 (0.2) 


0.9 (0.2) 


0.5 (0.2) 


0.0 (0.2 


0.41 (0.04) 


0.83 (0.02) 


0.56 (0.02) 


0.46 (0.02) 


this paper 



Re fer ences. — flUMaund et al.ll2007bl . (2) ITanaka et al.ll2009al (3) IMaund et al.ll2009l . (4) IKawabata et al.|[200l . (5) IWang et al~ll2003bl (6) ILeonard et al.l 
120021 . f7 HTanaka et al.H20o3 



